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We propose to perform a combined analysis of B — > tttt and B s — > K + K~ modes, in the framework 
of a global CKM fit. The method optimizes the constraining power of these decays and allows to 
derive constraints on NP contributions to penguin amplitudes or on the B 3 mixing phase. We 
illustrate these capabilities with a simplified analysis using the recent measurements by the LHCb 
Collaboration, neglecting correlations with other SM observables. 

CP violation in Bd, s decays plays a fundamental role in testing the consistency of the Cabibbo-Kobayashi-Maskawa 
(CKM) paradigm in the Standard Model (SM) and in probing virtual effects of heavy new particles. With the advent 
of the B-factories, the Gronau-London (GL) isospin analysis of B& — > tttt decays pQ has been a precious source of 
information on the phase of the CKM matrix. Although the method allows a full determination of the weak phase 
and of the relevant hadronic parameters, it suffers from discrete ambiguities that limit its constraining power. It is 
however possible to reduce the impact of discrete ambiguities by adding information on hadronic parameters [H [3] . 
In particular, as noted in refs. [1HS], the hadronic parameters entering the Bd — > tt + it~ and the B s — > K + K~ decays 
are connected by U-spin, so that the experimental knowledge of B s — > K + K~ can definitely improve the extraction of 
the CKM phase with the GL analysis. Indeed, in ref. [3], the measurement of BR(£> S — > K + K~) was used to obtain 
an upper bound on one of the hadronic parameters. 

After the pioneering studies performed at the TeVatron, very recently LHCb opened up the road to CP violation 
in B s — > KK decays [7j- The present experimental information is summarized in Table U] At present, one has all the 
necessary information to use the U-spin strategy proposed by Fleischer (F) in refs. [U [5] to extract the CKM phase 
from a combined analysis of Bd — > tt + it~ and the B s — > K + K~ decays. However, as we will show explicitly below, 
this strategy alone suffers from a sizable dependence on the breaking of U-spin symmetry [5] . 

Furthermore, in the B s system the measurement of any time-dependent CP asymmetry cannot be directly translated 

into a measurement of the angle /3 S = arg ^— pJ^J , even in the case of the so-called "gold-plated" b — > ccs decays. 

This is due to the fact that the angle (3 S is small and correlated to the subdominant amplitude in b — > ccs decays. 
Thus, measuring f3 s requires the determination of the subdominant decay amplitude. This is evident by noting that 
using CKM unitarity the b — > s decay amplitude can be written as 

A = V ub V: s T + V cb V*P or A = V u bV* s (T - P) + V t bV* s (-P). (1) 

Naively dropping the doubly Cabibbo-suppressed term proportional to V u bV* s would lead to the conclusion that the 
CP asymmetry measures 2f3 s with the first choice or that the CP asymmetry should vanish in the second choice. 
Clearly, a full treatment of the decay amplitude, taking into account correlations between the various CKM terms, is 
necessary to give a meaningful interpretation to the CP asymmetry. This is at variance with the Bd case, where the 
angle /3 is large and thus the time-dependent CP asymmetry in b — > ccs decays gives sin2/3 with a good accuracy^ 
In this respect, the combined analysis of the GL modes and B s — > K + K~ is optimal, since one has full knowledge 
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TABLE I. Experimental data used in the analysis. The correlation column refers to the S and -4cp measurements. Except 
from the results in ref. \7\, all other measurements have been averaged by HFAG [24]. The CP asymmetry of B + — > ty + tv° has 
been reported for completeness, although it has not been used in the analysis. 



1 The uncertainty due to the subleading amplitude can be quantified using Bd — > J/fp^, an d the correlation with the CKM terms in the 
decay amplitude is negligible l9lll2| 
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of the U-spin related control channel B d -> tt + 7t~, similarly to the case of B s -» fW°F"W° and _B d — s- i^(*) j^"(*) 
proposed in ref. [13]. Conversely, the "gold-plated" £? s — > J/ip4> decay has no U-spin related control channel, making 
the extraction of fi 8 problematic [Ti] . 

We propose to perform a combined analysis of the GL modes plus B s — > K + K~, including the time-dependent 
CP asymmetries, to obtain an optimal determination of the CKM phase within the SM. We show that this combined 
strategy has a mild dependence on the magnitude of U-spin breaking, allowing for a solid estimate of the theory error. 

Beyond the SM, NP can affect both the B d s — B d s amplitudes and the b — > d, s penguin amplitudes. Taking 
the phase of the mixing amplitudes from other measurements, for example from b — > ccs decays, one can obtain a 
constraint on NP in b — > s penguins. Alternatively, assuming no NP in the penguin amplitudes, one can obtain a 
constraint on NP in mixing. 

In this letter, we illustrate the points above in a simplified framework, neglecting SM correlations with other 
observables and using as input values sin2/3 = 0.679 ± 0.024 [21] and 2/3 s = (0 ± 5)° [25], obtained from b — > ccs 
decays. Clearly, the optimal strategy will be to include the combined analysis of the GL and F modes in a global fit 
of the CKM matrix plus possible NP contributions. 

The GL and F analyses were formulated with different parameterizations of the decay amplitudes. In order to 
use the constraints in a global fit, one should write the decay amplitudes with the full dependence on CKM matrix 
elements, but for the present analysis we can choose the F one and write the amplitudes as follows: 



A(B d -> 7T+7T-) = C(e*T - de ie ) , A(B d -> tt+tt-) = C(e~^ - de ie ) , 

A(B d -> vr 7r ) = ^=(Te l6T e^ + de l9 ) , A{B d -> 7r°7r°) = ^(Te^e^ + de l9 ) , 

A(B+ -> TT+TT ) = A(Bd ^ + "' ) + A(B d -► TT°TT°) , A(B~ -)■ T^TT ) = + A(B d -> TT TT°) 

A(B S -> K+K-) = C" T= A 75 (e^ + ±#dV e ') , A(B S -> F+F~) = C" I= ^ 75 ( e -^ + ±#dV e ') 



(2) 



where the magnitude of V u bV ud has been reabsorbed in C, the magnitude of V c bV cd / (V u bV ud ) has been reabsorbed in d 
and A = 0.2252. In the exact U-spin limit, one has C — C , d — d! and 8 — 8' . We have neglected isospin breaking in 
B d — > tttt, since its impact on the extraction of the weak phase is at the level of 1° [2"6H2l?] . The physical observables 
entering the analysis are: 

BR(B -> MM) = F{B) \MB^MMr + \A { B^MM) ? (3) 
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where F(B d ) = 1, F(B+) = T B+ /r Bd = 1.08, F(B.) = *(B.)(2 - (1 - y 2 )r(B s -> K+K-)/t b .) 00], t Bs = 



(1.425 ± 0.041)ps, $(B S ) = T Ss /T Bd (m 2 Jd /m| s ) v /(M 2 s -4M| r+ )/(M 2 d -4M 2 + ) = .9112, y s = Ar s /(2r s ) = 

(0.149 ± 0.015)/2 [31 , t(B s -> FJ+FJ") = (1.463 ± 0.042) ps [33 [33] and </> M (-B d ) = 2/3, <j) M {B s ) = -2/3 s in the 
SM. 

In the GL approach, one extracts the p.d.f. for the angle a = 7r — f3 — 7 of the Unitarity Triangle (UT) from the 
measurements of the three BR(B —> tttt), S(tt + it~), Acp{tt + it~ ") and ^ C p(7r 7r°)Q In this way, a (or, equivalently, 
7), is determined up to discrete ambiguities, that correspond however to different values of the hadronic parameters. 
As discussed in detail in ref. [3J, the shape of the p.d.f. obtained in a Bayesian analysis depends on the allowed 
range for the hadronic parameters. For example, using the data in Table [TJ solving for C and choosing flat a-priori 
distributions for d £ [0, 2], 8 € [— tt, 7t], T € [0, 1.5] and 9t G [— tt ; t] we obtain the p.d.f. for 7 in Fig.[l] corresponding 
to 7 = (68 ± 15)° (7 € [25, 87]° at 95% probability). Here and in the following we plot 7 only in the range [0, 180]° 
since the result is periodic with period 180°. 

Using instead the F method, one can obtain a p.d.f. for 7 from BR(B — > Tr + 7r~), BR(B S —> K + K~), S(ir + ir~), 
■^■Cp( 7r+7T ^)^ S{K + K~) and _4cp(F+ K~) given a range for the U-spin breaking effects. Fleischer suggested to 
parameterize the U-spin breaking in C'/C using the result one would obtain in factorization, namely 

rfact = \C'/C\ iact = 1.46 ±0.15, (4) 



2 Using unitarity of the CKM matrix, it is possible to write the B — > 7T7t decay amplitudes and observables in terms of a instead of 7 and 
0. However, for the purpose of connecting B — ► 7T7r to B s — > KK it is more convenient to use the parameterization in eq. 111. 
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FIG. 1. From left to right: P.d.f. for 7 obtained using the GL method as described in the text; p.d.f. for 7 obtained using the 
F method for k = 0.1, 0.3, 0.5. Here and in the following, dark (light) areas correspond to 68% (95%) probability regions. 




FIG. 2. From left to right: P.d.f. for 7 obtained using the combined method for n — 0.1, 0.3, 0.5; 68% probability region for 7 
obtained using the combined method (filled area) or the GL method (horizontal lines) as a function of k. 



where we have symmetrized the error obtained using light-cone QCD sum rule calculations in ref. |34j . However, this 
can only serve as a reference value, since there are nonfactorizable contributions to C and C that could affect this 
estimate [5]. In our analysis, we parametrize nonfactorizable U-spin breaking as follows: 

C = r lact r c C , d'e w ' = de w + r d de lr » , (5) 

with rc, rd and rg uniformly distributed in the range [1 — k, 1 + «], [0, k] and [— 7T, it] respectively. 

In Fig. [l] we present the p.d.f. for 7 obtained with the F method for three different values of k. We see that the 
method is very precise for k = 0.1, it is comparable to the GL method for k = 0.3, and it becomes definitely worse 
for k = 0.5. Thus, a determination of 7 from the F method alone is subject to the uncertainty on the size of U-spin 
breaking. 

We now consider the result of the combined GL+F analysis. In Fig. [2] we present the p.d.f. for 7 for k = 0.1, 0.3 
and 0.5. We see that the result of the combined analysis is much more stable against the amount of U-spin breaking 
allowed. We also plot the 68% probability region for 7 obtained using the combined method as a function of k, 
and compare it to the GL result. We see that there is a considerable gain in precision even for gigantic values of k. 
Actually, as can be seen in Fig. [3j where the posteriors for hadronic parameters and the U-spin breaking parameter 
rc are reported, the 68% probability range for rc is between ~ 0.4 and ~ 0.9. The fact that the rc posterior is not 
centered around 1, but the product rcr{ ac t is close to 1, may signal a failure of factorization and/or of the QCD sum 
rule estimate of rf act . On the other hand, the posteriors for d' and 9' are well compatible with small U-spin breaking. 
In any case, we think that the lesson to be learned from Fig.[3]is that values of k up to 0.6 or 0.7 cannot be excluded, 
but nevertheless the combined method remains useful. This happens because the peak around 7 ~ 30° in the GL 
result corresponds to values of 6 that are different from the ones needed in the F analysis to obtain similar values of 
7, while the peak at 7 ~ 70° is obtained for the same values of hadronic parameters in both the GL and F analyses. 

New Physics could affect the determination of 7 in the combined method by giving (electroweak) penguin contri- 
butions with a new CP-violating phase. Let us assume for concreteness that NP only enters b -> s decays, so that the 
isospin analysis of the GL channels is still valid. In the framework of a global fit, one can simultaneously determine 
7 and the NP contribution to b — > s penguins. For the purpose of illustration, we can just use as input the value of 
7 from tree-level processes, 7t r oo = (76 ± 9)° [ST\, and look at the posterior for 7 and for the NP penguin amplitude. 
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FIG. 3. From left to right and from top to bottom: P.d.f. for C, d, 8, rc, d', 9' obtained using the combined method for 
k = 0.9. 




FIG. 4. From left to right: P.d.f. for 7, 0npi ^np an d #np obtained using the combined method for k = 0.5. 



Taking for simplicity equal magnitudes for the NP contribution to A(B S — > K + K ) and A(B S — > K + K ), we write 

A(B S -> K+K') = C' 1 _ X x2/2 (e 11 + ^^(d'e 16 ' + e^d' NP e w ^)) , (6) 
A(B S K+K~) = C ' 1 _ X X 2 /2 ^ + ^-( d ' e%6 ' + e-^' NP d' NP e ie '-r)) . 

Taking uniformly distributed d' NP & [0,2] and <^np;^np [ — 71 ^ 7r ] we obtain the p.d.f. reported in Fig. [4] for k — 0.5. 
It yields 7 = (74 ± 6)°, and a 95% probability upper bound on around 1. The bound is actually much stronger 
for large values of (f)' NP . 

Finally, we notice that B s — > KK decays can also be used to obtain information on 4>m(B s ). The optimal choice in 
this respect is represented by B s -> if(*)°£(*)° (with B d -> K { *^K^° as U -spin related control channel to constrain 
subleading contributions), since in this channel there is no tree contribution proportional to e 17 [13j . However, the 
combined analysis described above, in the framework of a global SM fit, can serve for the same purpose. To illustrate 
this point, we perform the GL+F analysis not using the measurement of 2/3 s from b — > ccs decays. In this way, we 
obtain 2f3 s — (6 ±14)° for k — 0.5. With improved experimental accuracy, this determination will become competitive 
with the one from b — ¥ ccs decays, since the theoretical uncertainty can be estimated more reliably in the case of 
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B s — ^ K + K~ decays, waiting for time-dependent analyses of the — > K^°K^° channels^] To illustrate the 
potential of this method, we have repeated the analysis reducing the experimental uncertainty on Acp(Bd — > it + tt~), 
S(B d -> 7T+7T - ), A C p(B s K + K~) and 5(.B S -> K + K~) down to ±0.02. With such small experimental errors, it 
becomes crucial to take correctly into account the effect of the subleading term proportional to e* 7 in the amplitude 
(this is the case for any channel used to extract j3 s with an uncertainty of few degrees, including B s — > J/ip(f>). 
This is best done in the context of a global fit. For the purpose of illustration, we take as input the SM fit result 
7 = (69.7 ±3.1)° [31] and obtain 2f3 s — (2.6 ±2.7)° for k = 0.5. The error, which includes the theoretical uncertainty, 
could be further reduced improving the other relevant measurements, including the B^ decay modes, and by adding 
the Bd, s — > K^*'°K^*'° channels, allowing to test the SM prediction for CP violation in B s mixing. 

To conclude, let us summarize our findings. We suggest that the usual GL analysis to extract a from Bd — > tttt 
be supplemented with the inclusion of the B s — > K + K~ modes, in the framework of a global CKM fit. The method 
optimizes the constraining power of these decays and allows to derive constraints on NP contributions to penguin 
amplitudes or on the B s mixing phase. We have illustrated these capabilities with a simplified analysis, neglecting 
correlations with other SM observables. 
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